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^ '. Abstract 

We review the proposal of a quantum algorithm for Hilbert's tenth problem and provide further 
arguments towards the proof that: (i) the algorithm terminates after a finite time for any input 
of Diophantine equation; (ii) the final ground state which contains the answer for the Diophantine 
equation can be identified as the component state having better-than-even probability to be found 



by measurement at the end time-even though probability for the final ground state in a quantum 
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adiabatic process need not monotonically increase towards one in general. Presented finally are 



(T) . the reasons why our algorithm is outside the jurisdiction of no-go arguments previously employed 

o . 

to show that Hilbert's tenth problem is recursively non-computable. 
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I. INTRODUCTION 



We have claimed 0,0,0,0] to have a quantum algorithm for Hilbert's tenth problem 0] 
despite the fact that the problem has been proved to be recursively noncomputable. The 
algorithm makes essential use of the Quantum Adiabatic Theorem (QAT) p and other 



results (to be 
Computation 



urther substantiated in this paper) in the framework of Quantum Adiabatic 
?J in order to provide a single, universal procedure, physical or otherwise, 
which consumes finite amount of resources and terminates in a finite time and which can 
tell, in principle, whether any given Diophantine equation has any non-negative integer 
solution or not-thus solving Hilbert's tenth. 

We refer the readers to the above references and elsewhere for the history of Hilbert's tenth 
problem and its paramount importance in Mathematics and Theoretical Computer Science, 
as well as in Philosophy. In the next Section, we summarise the algorithm which employs 
Quantum Adiabatic Processes (QAP). Following that, we provide some arguments asserting 
why there is no level crossing, as required by the QAT, in the spectral flow associated 
with the QAP, except possibly at the end points. This will ensure that our algorithm can 
terminate in a finite time duration which is dictated by the smallest, non- vanishing energy 
gap separating the instantaneous ground state and relevant excited state in the spectral 
flow. We next provide the proof for the criterion identifying the final ground state of two- 
state systems, and then argue that this criterion is also generalisable to systems of infinitely 
many states of our algorithm. We also consider the algorithm in the context of existing 
no-go arguments for Hilbert's tenth problem in order to point out that the algorithm is 
outside the jurisdiction of those arguments. The paper is concluded with some remarks. 

II. THE QUANTUM ADIABATIC ALGORITHM 
A. Statement of the algorithm 

We first introduce the occupation-number state \n), n = 0, 1, 2, . . ., and the creation and 
annihilation operators and a respectively 

a|0) = 0, 

a\n) = y/n\n — 1), 
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o)\n) = \/n + \\n + \), 
(a'a)|n) = n\n), (1) 

and the coherent state, with complex number a, 

n=0 

Now, given a Diophantine equation with K unknowns, 

D(x 1 ,...,x K ) = 0, (3) 

we provide the following quantum algorithm to decide whether this equation has any non- 
negative integer solution or not: 

1. Construct /simulate a physical process in which a system initially starts with a direct 
product of K coherent states 

\m) = iw,) = (g)K>, (4) 
%=\ 

and in which the system is subject to a time-dependent Hamiltonian over the time 
interval [0, T], for some time T, 

with the initial Hamiltonian 

K 
i=l 

and the final Hamiltonian 

Hp = (D(a\a 1 ,...,a\ c a K ) y ) • (7) 

2. Measure/calculate (through the Schrodinger equation with the time-dependent Hamil- 
tonian above) the maximum probability to find the system in a particular occupation- 
number state at the chosen time T, 

P(T) = max|(^(T)|M)| 2 , 
IW> 

= \(m\{n}o)\\ (8) 
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where |{n}o) (which is a direct product of K particular occupation-number states, 
®f = i [n-i)) provides that maximal probability among all other direct products of K 
occupation-number states. 

3. If P(T) < 1/2, increase T and repeat all the steps above. 

4. If P(T) > 1/2, then |{n}o) is the ground state of Hp (assuming no degeneracy-we 
discuss the case of degeneracy below) and we can terminate the algorithm and deduce 
a conclusion from the fact that: 

H P \{n} ) = iff the equation (J3J) has a non-negative integer solution. (9) 



We refer to 



□ □ □ 



for the motivations and discussions leading to the algorithm above. 



Note that it is crucial that Hj does not commute with Hp, 

[Hi, Hp] + 0, (10) 

and of course both Hamiltonians are of infinite dimensions. (This condition can be relaxed 
to that of unbounded dimensions, see later.) 

To remove any possible degeneracy of the ground state for any Hp, we can always in- 
troduce to Hp a symmetry-breaking term of the form (7a] + 7* a*) for some % in the limit 
I7I — ► |4|. This term destroys the symmetry generated from the commutation between 
Hp and the occupation-number operators. However, we can always recover the symmetry 
in the limit and modify the algorithm above slightly to reach an answer for the Diophantine 
equation in question. Note also that the arguments of the Section IIHI below can be easily 
modified to accommodate this symmetry-breaking term to ensure there is no ground-state 
degeneracy inO<t/T<l. 



B. Its probabilistic nature 

We will prove in the next Section that there exists a finite T whence the algorithm can 
be terminated. But it is important to recognise hereby that our algorithm is probabilistic 
in nature. As the halting criterion, the probability P{T) of JHJ) for certain state need to be 
more than one-half. If the algorithm is implemented by certain physical process then we can 
approximate this probability by the relative frequency that a particular state is obtained in 
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repeating the measurements many times over. The Weak Law of Large Numbers, however, 
can only assert that this frequency is within a distance e of P{T) with certain probability of 
at least 1 — 5. It is this probability 1 — 5 that gives the quantum algorithm its probabilistic 
nature. Both e and 5 are dependent on the number of measurement repetitions from which 
the frequency is obtained-see, for example, |8| for a statement and proof of the theorem. In 
general, we can always reduce e and 5 to arbitrarily close to zero by increasing the number 
of repetitions L appropriately, 

L > l/(4e 2 5). (11) 

As an anticipation of the numerical simulations of the algorithm on recursive computers 0, 
, we mention now that the probabilistic nature of the algorithm is manifest differently 
there through the extrapolation to zero step size in solving the Schrodinger equation. Such 
extrapolation is necessary for a correct truncation, for a given Diophantine equation, of the 
dimensions of the Hilbert space involved. 

We will later come back to this probabilistic nature in Section IS in the discussion of 
no-go arguments for Hilbert 's tenth problem. In ending this Section, we wish to point out 
here that, contrary to common misunderstanding, probabilistic computation in general is 
not equivalent in terms of computability to Turing computation 



III. NO LEVEL CROSSING FOR THE GROUND STATE 

One of the sufficient conditions for the Quantum Adiabatic Theorem 

[6] is that the 

relevant level, which in our case is the instantaneous ground state, doesnot cross with any 



other level as time progresses. In 
the time-dependent Hamiltonian 



4] 



we have given some arguments for no level crossing for 



H{sT) = (l-s)H I + sH P , (12) 

with the reduced time s in the interval s e (0, 1). Now we present below some other ar- 
guments adapted from those of Ruskai 111, which in turn are based on Perron- Frobenius's 

n 

theorem |12| for finite but unbounded number of dimensions, and on Reed-Simon's theo- 

n 

rems [13J for infinite number of dimensions. 
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Let us consider the exponential operator exp(— TC(sT)), it can be expressed by the Lie- 
Trotter product formula as 



-H(sT) 



lim fe-sH^-TT Y,M a i+< a i) 



(13) 



For general operators in dimensionally infinite Hilbert spaces, the convergence in the above 



is understood to be the strong convergence of the rhs to the lhs 



2l| 



A. Finite but unbounded number of dimensions 



Our algorithm of the last Section actually requires only sufficiently large but finite 
Hilbert space with a truncated basis consisting of occupation-number vectors, {\n) : n = 
0, 1, ... , N}, for some N beyond which higher occupation- number states contribute negligi- 
bly to the dynamics of low-lying states relevant to our problem. This follows from the facts 
that the normalised wavefunction at any given time has a support which spreads signifi- 
cantly only over a finite range of occupation-number states, and that the explicit coupling, 
and thus the influence, between one particular instantaneous eigenstate and other dimin- 
ishes significantly outside some finite range of occupation-number states (as can be seen 
through a set of differential equations in |4| connecting instantaneous eigenstates oi"H(t) at 
different time t). We will exploit these facts fully in the numerical simulations of the quan- 
tum algorithm . How large is sufficient, and thus how much the truncation N should 
be, depend very much on the particular Diophantine equation being investigated. Such 
truncation will be discussed and implemented in the numerical ssimulations of the quantum 
algorithm 

a a. 

In the last Section, we need not worry about the size of N as we have 
employed dimensionally infinite Hilbert space. However, for any arbitrary and unbounded 
N, the Hamiltonian (JT2Jl above has a representation which is a finite square matrix and in 
which all the off-diagonal elements are contained in Hj. 

The elements of the matrix product in the square brackets on the rhs of Q13jl . can be 
expressed in the truncation to an arbitrary N (where a^\N) =0) as 

JV 

(m|e fta *V* a< |n) = ^(m|e fta * |g)(g|e^ a! |n>, 



q=0 

N oo 

q=0 k=0 



n), 



(14) 



1=0 



where 

A = (15) 

It then follows that for {ctj} real and positive, and for < s < 1, the lhs above is always 
non-zero and positive. For example, when N > m > n, at least the term with k = m — n 
and I = will contribute nonvanishingly and positively to the sum of non-negative terms on 
the rhs to ensure the matrix element on the lhs is always positive. 

In the basis of occupation-number states, the first factor in the bracket of the rhs of (jlHj) 
is a diagonal matrix, whose net effect is to multiply each row of the matrix product in 
square brackets in (jl3j) by non-zero, positive numbers (being exponentials). As a result, for 
Q!j = a* > 0, for i = 1, 2, . . . , K, the product in the bracket on the rhs of (fE?j) has only 
positive elements, and so has its m-th power. Thus, Perron-Frobenius's theorem Q| for 
(finite) matrices with strictly positive elements can be applied to the lhs of (|T3| to confirm 
that the largest eigenvector of e~ n( - sT ^ is unique for these values of {c^}. That is, H(sT) 
itself has a unique ground state for < s < 1 (since, by construction, 7^(0) = Hi has a 
nondegenerate ground state). 

At s = 1, the matrix e _i ^r Hl becomes the identity matrix which has zero non-diagonal 
elements and becomes reducible, violating the conditions of Perron-Frobenius's theorem and 
thus allowing the possibility that 7i(T) = Hp has a degenerate ground state. 



B. Infinite number of dimensions 



To handle the infinite dimensions of the operator product in the square brackets of (|TH|) 
directly, that is, without any truncation, we can employ the so-called holomorphic represen- 



tation 
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3 

m\( 1 • < 1 "'|n) 



to arrive at 



/W /3?aii_\ / d V d V dujduj d£d£ 



27ri 2iri 2iri ym\ 

i fd\ m ( dXU 

VmM \ds) \ds 

i fd\ m fd xn 



drjdr] dC, d^ 
2ni 2m 



(16) 



s=s=0 



where we have integrated over {u, u} and introduced the source terms s and s in the second 
line, and done a Gaussian intergration to get the last line. It can be seen from (fTBj) . that the 



matrix element always is strictly positive for any positive integers m and n provided that 
< s < 1 and on is real and positive. From this we are led to the conclusion that 



[n e M 



g M S ^i aia iQ M Ei a i a i 



\m) 



EiOn+KI )^ n | e -jrE<«iOi e -5rEia*oi| m ^ > o, (17) 

and consequently that the lhs of (fT3*j) has positive matrix elements for real and positive ol{. 

For the case of bounded operator on the lhs of (j!3jl . we can adopt a theorem of Reed and 
Simon 2^, a cut-down version of which can be rephrased for our purpose as 

7i has a nondegenerate ground state iff e~ aH is positive improving for all a > 0. 

In our occupation-number basis, having e~ aH positive improving for all a > is equivalent to 
having matrix element (m\e~ an \n) strictly positive for all a > and for all positive integers 
m, n. But the nonvanishing of such matrix elements follows simply and directly from the 
result in (|T7j) and the strong convergence in (jl3)l . Thus, finally, the non- degeneracy of the 
ground state of the dimensionally infinite operator 7i(sT) also follows. 

We can in fact prove a stronger result, by invoking another theorem of Reed and Si- 
mon {23] ]. that all the eigenvectors of H(sT), not just the ground state, are non-degenerate 
for real and positive a^s and for < s < 1. 

C. Extension to a>i G C 

Having established the above results for real and positive a^s, we now extend them to 
all complex- valued but nonvanishing a^'s. Indeed, having all a« vanishing would make Hj 
commute with Hp, leading in general to an unwanted crossing of the ground state of 7~C(t) 
at some t G (0, T). 

The crucial point is to note that the Hamiltonian (jHJ) is invariant under the following 
phase transformations, for real and arbitrary 6t, 



if) 



e i6i ai. 



(18) 



Denoting 

a, = e ie ' ai = (V ie < a M a, [>* a M , (19) 



we see that these operators also satisfy the canonical commutation relation, [a^, aj] = 5ij, and 
we can thus construct a unitarily related Fock space, with the basis consisting of |{n<i}) = 
<S>f =1 \n.i), such that 

\fk) = e- in ' e '\ ni ), (20) 
(a\ai)\ni) = n^fii). 

Now, with nonvanishing complex- valued = e _l6i |ai|, we can perform the phase trans- 
formation (jl8j) to obtain an equivalent Hamiltonian containing only real and positive |aj|, 
and a%. The arguments of the two subsections above, for both finite and infinite number 
of dimensions, will carry through but this time with the unitarily transformed occupation- 
number states Once again, the nondegeneracy of the ground state of Tt(t) is assured, 
but this time extended for nonvanishing complex- valued a^s. 

All of the above agrees with our alternative arguments in [4] which lead to the nondegen- 
eracy of the ground state of TC(sT) for < s < 1. The possible degeneracy of Hp can then 
be handled as discussed in the last Section. It follows next from the Quantum Adiabatic 
Theorem that if the instantaneous ground state never crosses with any other state for the 
whole of the spectral flow then it takes a non- vanishing rate of change 1/T, and thus only 
a finite time T, for the probability of the final ground state to evolve arbitrarily closed to 
one. This implies that our quantum algorithm can always be terminated in a finite time. 
The next Section gives the criterion for finding that terminating time. 

IV. IDENTIFYING THE GROUNDSTATES 
A. The overall picture 

The crucial step of any quantum adiabatic algorithm is the identification of the ground 
state of the final Hamiltonian, Hp. Normally it is identified as the probabilistically dominant 
state obtained for an adiabatic evolution time T sufficiently long as asserted by the Quantum 
Adiabatic Theorem jfj. In our case we do not in advance know in general how long is 
sufficiently long (the Theorem offers no direct help here); all we can confidently know is that 
for each Diophantine equation and each suitable {a{\ there is a finite evolution time (the 
finiteness is due to the non- degeneracy of the instantaneous ground state-see last Section) 



9 



after which the adiabatic condition is satisfied. We thus have to find another criterion to 
identify the ground state. 

We will make full use of the fact that the eigenstates of Hp are by construction just the 
occupation-number states, = ® i=1 \rii) : (a|aj|rij) = 7ij|nj})&(ni = 0, 1, . . .)}, among 

which is the final ground state to be identified. 

The identification criterion we have found can be stated as: 

The ground state of Hp is the component state |{^}o) whose measurement 
probability is more than ~ after the evolution for some time T of the initial 
ground state \{a}j) according to the Hamiltonian ilfy) . 

\(ip(T)\{n}o)\ 2 > -, for some T, if and only if |{n}o) is the ground state of Hp. (21) 

We can recognise the Quantum Adiabatic Theorem in the 'if part of the statement above; 
and we need only prove the 'only if part. In general, it suffices that the initial ground 
state \gi) of some Hi should not have any dominant component in the occupation-number 
eigenstates \{n}) of H P , 

l<SilW>i 2 <^VW. (22) 

Our choice of initial coherent state (JU clearly satisfies this condition. 

Consequently, we only need to increase the evolution time T until one of the occupation- 
number states is obtained at time T with a probability of more than ~ then this will be our 
much desired ground state. 

We will prove this criterion for two-state systems in the next Subsection and then argue 
that it is also applicable for systems of finitely many and of infinite number of states because 
among those states only two states, which may or may not be the instantaneous ground state 
and first excited state, become dominantly relevant at any instant of time- provided we start 
out with the ground state of the initial Hamiltonian Hi. 

B. Two-state systems 

We divide the proof for two-state systems in turn into three parts: 
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• In the first one, we show (Eq. (|39|) below) that the maximum probability for the state 
which is not the final ground state cannot be more than | for any evolution time T- 
subject to certain condition stated below (and which is also satisfied by our choice of 
Hi in the generalisation to dimensionally infinite Hilbert spaces). 

Note that the probability need not be monotonic function of T. 

• In the second part, we appeal to the Quantum Adiabatic Theorem (f| which asserts 
that eventually for sufficiently large T (i.e. for sufficiently slow evolution rate 1/T) 
the probability for the final ground state approaches one. 

• Then by combining the two parts above, we can, without the need of knowing a priori 
how slow is sufficiently slow for the evolution rate 1/T, conclude that eventually the 
probability of one state will rise above | as T increases and that state must be the 
ground state. 

It suffices to establish the first part of the arguments above. Let \g(t)) and |e(t)) be 
the instantaneous eigenstates of W 2 '(t), the two-dimensional counterpart of (J3J). They are 
related to the eigenstates of HS 2 \t + St) by 

\g{t)) = cos((3(t)/2)\g(t + 6t))+sm((3(t)/2)\e(t + 6t)), 

\e(t)) = - sin ((3 (t) / 2) \g (t + 5t)) + cos ((3 (t)/ 2) \e(t + 6t)), (23) 



where 



tan(3(()) = mmm> - iPHf' + 0{H ' (24) 



by recalling that H {2) (t + 6t) = K {2 \t) + 5td t K {2) + 0{5t 2 ), where 



l p ) 



Hf ] = e 9 |^(0))(^(0)| + e e |e(0))(e(0)|, (25) 
Hf = r g \g(T))(g(T)\ + r e \e(T))(e(T)l 

for some e's and T's and such that (fTUjl is observed, 

[Hf\Hf] ^ 0. (26) 
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Note that by redefining the phases of \e(t)) and \g{t)) in (j23J) we can always choose 

< p(t)/2 < tt/2. (27) 

Assuming that at time (t + St), 

m + St)) ~ cos(0(t))|^(t))+e-^sin(0(t))|e(t)), (28) 
for some 0(t) and such that 

tan(0(t)) > 0, (29) 

then 

\iP(t + 28t)) = e- in{t+st)st \^(t + 6t)}, 

~ (cos(0(t)) cos(/?(t)/2) - e"*® sin(0(t)) sin(/5(t)/2)) \g(t + St)) (30) 
+e _ ie(4+(5t ) ( cos (0( t)) sin(/?(t)/2) + e- i9W sin(0(t)) cos(/?(t)/2)) |e(t + 5t)), 

by the use of (|2*%|) and (j2Hl) - According to (jHOJ), the probability to be found in the ground 
state at time (t + 2St) is 

cos 2 (0(t + St)) = |cos(0(t)) cos(/3(t)/2) - e~ im sin(0(t)) sin(/?(t)/2) | 2 , 

> cos 2 (0(t)+/3(t)/2); (31) 

where we have made use of (|2*nj) to derive the last inequality. Similarly, the probability to 
be found in the excited state is 

sin 2 (0(t + St)) = |cos(0(t)) sin(/3(*)/2) + e~ im sin(0(t)) cos(/3(t)/2) f , 

< sin 2 (0(t)+/?(t)/2), 

/ i/« \ 

< sin 2 1^/3(^/2 J , (32) 

where we have used the arguments of induction and the infinitesimality in St of (3{n5t) to 
arrive at the last line above j^j]. 
Denoting 

r/8t 

W ( T ) = J™ (33) 

n=0 



12 



we now sum up the infinitesimals 0(St) of (|24|) to have 

tmMT)) = L mmm - 5pw* (34) 

Our next step is to prove that the integrand in (JM)) is never zero for t G (0, T). It is a 
proof by contradiction if the opposite is assumed. That is, assuming that there exists a time 
to at which the numerator of the integrand vanishes 

= (e(t )\d t H^(t )\g(t )) = (e{t,)\Hf - tfj 2) |s(t ))/T, (35) 

by the use of (J25|) . This means that the combined operator Hp — Hj is diagonal in the 
basis {\g(t )), |e(t ))}. On the other hand, the operator ft (2) (t ) = Hf ] + (t /T)(H P 2) -Hf) 
is automatically diagonal in the same basis by virtue of its very definition. Thus, these two 
operators must commute 

= [H^{t Q ),Hf -Hf\ 

= [Hf\Hf], (36) 

which contradicts our assumption concerning Hi and Hp in lj'26j) . Thus there cannot exist 
any such t that satisfies (j3SJ)- 

Having established that the integrand is never zero, we can find out its sign by evaluating 
it at t = 0, say. The integral (|3~4"j) likewise never vanishes for any r and assumes the same 
sign as its integrand. 

The case of interest is when the sign is positive 

tan(ct>(r)) > 0, for all r G (0,T), (37) 

which implies that 

< uj(t) < 7r/2. (38) 

This last line follows from the continuity in r which smoothly connects w(r) at r ^ to 
w(0): 

• from (|2~7jl we know that /3(0) is in the range [0, 7r); 

• thus, uj(Q), which equals /3(0), has to be less than 7r/2 for (|3*7j) to hold at r = 0; 
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• consequently, as tan(u;(r)) never for any r is negative, uj(t) has to be less than tt/2, by 
the continuity of oo(t) as a function of r, for (pi?)) to hold at all r. (This last statement 
can also be proved by employing the method of proof by contradiction.) 

Thus, we are led from (jSHJ) and to an upper bound on the probability of the excited 
state at all time 

sin 2 (0(r)) < sin 2 (w(r)/2) < sin 2 (7r/4) = -. (39) 

This completes our main arguments for this Section. 

Note that the condition (|37|) is crucial here. It is satisfied for 

Ks(0)|e(T)>| 2 <l (40) 

and the system is initially in the ground state \i[>(0)) = Ifi'(O)). If the opposite of the above 
is true then we can show that there exists some range of T such that (f3H|) does not hold. 
In fact, this can easily be seen in the sudden approximation when the rate of change 1/T is 
large 

hm|^(T)> « hKo)) = to(o)>, 

from which, assuming the opposite of PU|). 

|WT)|e(T))| 2 « |(^(0)|e(T))| 2 >i (41) 

which is the opposite of the result (}39|) . (The existence of a range of values of such T 
immediately follows from the continuity of \ip(T)) as a function of T.) By solving the 
Schrodinger equations with time-dependent Hamiltonians, we illustrate in Figure 1 both the 
cases when (l3~9"j) is and is not held depending on whether the condition (jjOj) is satisfied or 
not. 

When the above is generalised to dimensionally infinite systems of the quantum algo- 
rithm for Hilbert's tenth problem, the counterpart of ()23)1 is the set of differential equations 
connecting the sets of instantaneous eigenstates of Ti.it) at different time t ^|. There, our 
choice of the coherent state (jlj), as the ground state in which the system initially has to be, 
entails that the condition (}4T)j) is always satisfied, since for any n and a 

,o\n>\ 2n 1 

|(a|n)| 2 = e"l Q l 4p<-. (42) 
14 




FIG. 1: Probability distributions as functions of T for two-state systems, illustrating the cases of 
(A) |(5(0)|e(T))| 2 = 3/4 > 1/2 and (B) \(g(0)\e(T)}\ 2 = 1/2. For the first case, there is a range 
of T which violates (|39|) . while in the latter we always have (j39[) . The excited-state probability 
always approaches zero for sufficiently large T according to the Quantum Adiabatic Theorem. 

And we expect that, in these infinite systems too, the probability to be subsequently found 
in any particular excited state cannot be greater than | at anytime. Such arguments have 
also been numerically confirmed in several simulations of our algorithm and also of a 
modified version 1 1 611 - 



V. HOW CAN WE COMPUTE THE NON-COMPUTABLE? 

When proposed the tenth problem in 1900, Hilbert himself never anticipated the link it 
would have with what is the Turing halting problem of the yet-to-be-born field of Theoretical 
Computer Science. The Turing halting problem was only introduced and solved in 1937 by 
Turing, and the link of equivalence between the two problems was only established in 1972 p. 



A. The Turing halting problem 

The question of the Turing halting problem can be phrased as whether there exists a 
universal process according to which it can be determined by a finite number of operations 
if any given Turing machine would eventually halt (in finite time) starting with some spe- 
cific input. Turing raised this problem in parallel similarity to the Godel's Incompleteness 
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Theorem and settled it with the result that there exists no such recursive universal proce- 
dure. The proof is based on Cantor's diagonal arguments, also employed in the proof of the 
Incompleteness Theorem. 

The proof is by contradiction starting with the assumption that there exists a recursive 
(and hence Turing computable) single-valued halting function h(p, i) which accepts two 
integer inputs: p, the Godel encoded integer for the Turing machine in consideration, and 
i, the Godel encoded integer for the input for p, 

{0 if p halts on input %\ 
(43) 
1 otherwise. 

One can then construct a program Turing (n) having one integer argument n in such a way 
that it calls the function h(n,n) as a subroutine and then halts if and only if h(n,n) = 1. 
In some made-up language: 

Program Turing 
input n 
10 call h(n, n) 

(44) 

if h(n, n) = goto 10 

stop 

end 

Let t be the Godel encoded integer for Turing; we now apply the assumed halting function 
h to t and n, then clearly: 

h(t, n) = if and only if Turing halts on n 

if and only if h(n,n) = 1, (45) 

from which a contradiction is clearly manifest once we choose n = t. 

The elegant proof above was only intended by Turing for the non-existence of a recur- 
sive halting function. Unfortunately, some has used this kind of arguments to argue that 
there cannot exist any halting function in general! We have pointed out elsewhere the 
fallacies in such use, and considered carefully the implicit assumptions of Cantor's diagonal 
arguments. See also the Subsection IV CI below. 
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B. The equivalence between Hilbert's tenth and the Turing halting problem 



It is easy to see that if one can solve the Turing halting problem one can then solve 
Hilbert's tenth problem. This is accomplished by constructing a simple program that sys- 
tematically searches for the zeros of a given Diophantine equation by going through the 
non-negative integers one by one and stops as soon as a solution is found. The Turing 
halting function (existed by assumption) can then be applied to that program to see if it 
ever halts or not. It halts if and only if the Diophantine equation has a non-negative integer 
solution. 

Proving the relationship in the opposite direction, namely that if Hilbert's tenth problem 
can be solved then will be Turing halting problem, is much harder and requires the so-called 
Davis-Putnam- Robinson-Matiyasevich (DPRM) Theorem j^|: 



Every recursively enumerable (r.e.) set r L 25] of n-tuple of non-negative inte- 
gers has a Diophantine representation. That is, for every such r.e. set there 
is a unique family of Diophantine equations D(ai, . . . , a n ; x%, . . . , x m ) = 0, each 
of which has n non-negative integral parameters (ax, ... , a n ) and some m vari- 
ables (x\, . . . , x m ), in such a way that a particular n-tuple (ax, ... , a n ) belongs 
to the set if and only if the Diophantine equation corresponding to the same n 
parameters has some integer solutions 

Let us number all Turing machines (that is, programs in some fixed programming language) 
uniquely in some lexicographical order, say. The set of all non-negative integer numbers 
corresponding to all Turing machines that will halt when started from the blank tape is 
clearly a r.e. set. Let us call this set the halting set, and thanks to the DPRM Theorem 
above we know that corresponding to this set there is a family of one-parameter Diophantine 
equations. If Hilbert's tenth problem were recursively soluble, that is, were there a recursive 
method to decide if any given Diophantine equation has any solution then we could have 
recursively decided if any Turing machine would halt when started from the blank tape. 
We just need to find the number representing that Turing machine and then decide if the 
relevant Diophantine equation having the parameter corresponding to this number has any 
solution or not. It has a solution if and only if the Turing machine halts. 

But that would have contradicted the Cantor's diagonal arguments for the Turing halting 
problem! Thus, one comes to the conclusion that there is no single recursive method for de- 




17 



ciding Hilbert's tenth problem. For the existence or lack of solutions of different Diophantine 
equations one may need different (recursive) methods anew each time. 

Now, having reached this conclusion, we wish to point out that logically there is nothing 
wrong if there exist non-recursive or non-deterministic or probabilistic methods for deciding 
Hilbert's tenth. 



C. The quantum algorithm in context 

In claiming that our quantum algorithm can somehow compute the noncomputable, we 
also need to consider it in the context of the no-go arguments above. Those arguments, 
indeed, cannot be applicable here because of several reasons. First of all, the proof for the 
working of our algorithm is not quite constructive, implying its non-recursiveness in some 
sense. The mathematical proof for the criterion for ground-state identification, as can be 
seen from the Section IIV[ is highly non-constructive as we have had to employ at several 
places the methods of proof by contradiction and also of (non-constructive) analysis for 
continuous functions. 

Secondly, and more explicitly, our algorithm is outside the jurisdiction of those no-go 



arguments because of its probabilistic in nature. We have argued elsewhere 10] against the 
common misunderstanding that probabilistic computation is equivalent in terms of com- 
putability to Turing recursive computation. They are not! We have pointed out that if a 
non-recursively biased coin is used as an oracle for a computation, the computation carries 
more computability than Turing computation in general. Here, we shall show explicitly in 



the below how the probabilistic nature of our a 



gorithm can avoid the Cantor's diagonal 



no-go arguments presented in Subsection IV Al [27 1. 

Because our algorithm is probabilistic (see Subsection 111 B|) . we can only obtain an answer 
with certain probability to be the correct answer. This probability can be, with more and 
more work done, made arbitrarily closed, but never equal, to one as shown in (jllj) . Thus, 
instead of the halting function (j43|) our algorithm can only yield a probabilistic halting 
function ph, similar but not quite the same as it must have three arguments instead of two, 

if p halts on input i, with maximun error-probability 5; 
ph(p,i,5) = { ' (46) 

1 if does not halt, with maximun error-probability 5. 
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In order to follow the Cantor's arguments as closely as possible, we need to restrict to 

5 = 2" J , (47) 

for some integer J. Following the flow of Cantor's arguments, one can then construct a 
program pTuring(n, S) having arguments n and S in such a way that it calls the function 
ph(n,n,S) as a subroutine and then halts if and only if ph(n,n,5) = 1. In some made-up 
language: 

Program pTuring 
input (n,6) 
10 call ph(n, n, 5) 

(48) 

if ph(n, n, S) — goto 10 

stop 

end 

Similarly, let t p be the Godel encoded integer for pTuring (barring the case pTuring does 
not have an integer encoding, which is quite possible for a quantum algorithm Q]). We 
next apply the probabilistic halting function ph to t p and (n, 5), which is the total input for 
pTuring, and with some 5', to obtain: 

ph(t p , h, 5') = iff pTuring halts on n, with maximum error-probability 5', (49) 

where h encodes (n, 5) uniquely, that is, h = Pipi, where p\ and P2 are two different prime 
numbers. No matter what we choose for n and 5' and J, we cannot diagonalise the above, 
because n ^ n, unlike previously (j45jl . Thus we never run into mathematical contradiction 
here. 

The above arguments apparently have nothing to do with being quantum mechanical, but 
everything to do with being probabilistic. Quantum mechanics, however, has given us an 
inspiration to realise a probabilistic algorithm capable of deciding Hilbert's tenth problem 
with a single, universal procedure for any input of Diophantine equations. 

VI. CONCLUDING REMARKS 

In this paper we have provided the analytical results and arguments for the working 
of our quantum algorithm for Hilbert's tenth problem. Numerical simulations for some 
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simple Diophantine equations have been preliminarily reported in |3[ and will be available 
fully elsewhere where we shall explain how to cope with the required dimensionally- 
unbounded Hilbert space and also argue that while the algorithm has been inspired by 
(quantum mechanical) physical processes it may be simulated on Turing computers, despite 
the proof that Hilbert 's tenth problem is recursively noncomputable. 
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